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Abstract 

In this paper, we propose two new models in f(T) gravity to realize universe acceleration and 
phantom crossing due to dark torsion in the formalism. The model parameters are constrained 
and the observational test are discussed. The best fit results favors an accelerating universe with 
possible phantom crossing in the near past or future followed respectively by matter and radiation 
dominated era. 
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1. INTRODUCTION 



In recent cosmological observations such as the Type la Supernova (SNe la) PQ[2], the 
cosmic microwave background radiation [3] [4] and the large scale structure jSjfB], et al., we 
found that in near past, the universe is undergoing an accelerating expansion phase. This 
surprising behavior represents one of the most complex issues in today cosmology. So far, the 
so called "dark energy" (DE) in our universe, as an exotic energy component with negative 
pressure dominating the universe is known as the most prominent candidate driving universe 
to an accelerating expansion phase. 

The simple ACDM model with a constant equation of state (EoS) is known as the 
simplest model in general agreement with the current experimental data, However, the model 
suffers from "cosmological constant problem [7] [8j. In addition, recent observation reveals 
that the model may endure from an age problem [9J. An alternative approach thus seems 
reasonable to find and investigate for the current acceleration of the universe. Observational 
probes support a small deviations in the EoS parameter by crossing phantum divide line 
from above to below in the near past. This therefore requires a dynamical description 
of the parameter [ID] . The accelerating universe may be interpreted as a failure in our 
understanding to the gravitational law. Thus we may require an alternative to standard 
theory of gravity. 

In last few years, a variety of cosmological models, such as, quintessence [IT] , k-essence 
[12], tachyon [13] [HJ, phantom [13], quintom [ID] [IB]; Chaplygin gas and its generalization 
[IT] [IE] , holographic DE P3]PDJ, agegraphic DE [22] [23], Ricci DE [2j] have been investi- 
gated. 

In addition, a modification in gravity such as in f(R) [25]-[35]. or other curvature invari- 
ants [36] [39], by coupling a scalar field to curvature [40J [38J [41J, a vector field contribution 
[42j . or in higher dimensional spacetimes, have widely been studied [43]. Some of these 
models that display both universe acceleration and phantom crossing are supported by ob- 
servational probes [44]— [48]. 

In a remarkably alternative approach to standard general relativity, we use Weitzenbck 
connection instead of Levi-Civita connection in the theory where subsequently curvature is 
replaced by torsion. In this approach the torsion is formed completely from products of 
only first derivatives of the tetrad. This theory known as " Teleparallelism" , was originally 
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introduced by Einstein in 1928 [49j[50j. It differs from standard theory of relativity only in 
"boundary terms" involving total derivatives in the action. The " Teleparallelism" theory is 
formulated by gauging external spacetime translation. It employs the Weitzenbck spacetime 
which is characterized by the metricity condition and vanishing of the curvature tensor. The 
approach contains a variety of distinctive manifestation both from physical and geometrical 
aspects [5T]~[65]. 

There are models based on the teleparallel equivalent of general relativity (TEGR) that 
presented as an alternative to inflationary models with no need to inflaton field [58j[59j. 
Also there are DE models constructed from teleparallelism that explain current acceleration 
of the universe without using DE component j6()]-[66] In these models dark torsion (DT) 
is responsible for the acceleration of the universe, and the field equations are always 2nd 
order equations. Among other advantages, this characteristic makes these theories simpler 
than the resulting in f(R) theories. Note also that in generalized theories of TEGR, where 
the Lagrangians are algebraic functions of the usual teleparallel Lagrangian, both actions 
and field equations are not invariant under local Lorentz transformations [2D]. However, the 
usual teleparallel theory equivalent to general relativity, is just a special case of the so-called 
f(T) gravity where the Lagrangian of teleparallel gravity, is substituted by T The f(T) 
gravity provides an alternative to conventional dark energy in general relativistic cosmology 
to explain the acceleration of the universe. 

Note that the main requirement of teleparallelism and f(T) gravity is that there exist 
a class of frames where the spin connection vanishes, but torsion does not [2T]. In fact, 
the torsion tensor is formed solely from products of first derivatives of the tetrad and the 
spacetime torsion manifests itself as a generator of gravitational repulsion. This enables the 
theory to produce universe acceleration (for example, see [67] [75j ) 

The original motivation for considering f(T) gravity is to account for the late time cosmic 
acceleration without the need for a cosmological constant and/or dark energy. Therefore, 
it makes sense to restrict ourselves to specific models that exhibit this property. In this 
manuscript, we follow the authors in [60] and introduce two new forms of f(T) models to 
investigate the universe acceleration and phantom crossing by best fitting the models with 
the observational data. 
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2. GENERAL CONSIDERATIONS 



Teleparallelism uses a vierbein field ej(x M ), i = 0, 1, 2, 3, as a dynamical object which is an 
orthonormal basis for the tangent space at each point x 11 of the manifold: e^e^ = rjij , where 
rjij = diag(l, — 1, —1, —1). The vector can be described by its components ef , /x = 0, 1, 2, 3 
in a coordinate basis; i.e. = ef <9 M . Note that Latin indices refer to the tangent space, while 
Greek indices label coordinates on the manifold. The metric tensor can be obtained from 
the dual vierbein as g^x) = r]ije l Jx)e{,{x) . Though general relativity uses the torsionless 
Levi-Civita connection, in Teleparallelism we apply the torsionless Levi-Civita connection, 
Teleparallelism uses the curvatureless Weitzenbock connection, whose non-null torsion is 

= f*. - 1%, = ef(0„4 - d v e\). (1) 

The above equation contains all the information about the gravitational field. The TEGR 
Lagrangian is constructed from torsion equation, 0, and the dynamical equations for vier- 
bein are corresponding to the Einstein field equations for the metric 
The teleparallel Lagrangian is 

T = S/TV, (2) 

where 

= \{K» V P + 5»T 6v e - 6;t\) (3) 
and K^ u p is the contorsion tensor 

K»% = ^(T% - T v % - T/"), (4) 

that is equal to the difference between Weitzenbock and Levi-Civita connections. 

In this work we start with the following action where Lagrangian density sums over T 
and f(T). Thus the action reads 

where e = det{e l ) = \f—g. The action only with T corresponds to TEGR. In the presence 
of matter field coupled to the metric, variation of action with respect to the vierbein yields 
the following equations [63] 

e-%(eSr)(l + fiT)) - e^^S^il + f'(T)) 

+Srd,(T)f"(T) + 1<(T + f(T)) = 47rGe/T/, (6) 
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where prime denotes differentiation with respect to T, = efSp^ and T^ v is the matter 
energy-momentum tensor. 



3. COSMOLOGICAL SOLUTION AND OBSERVATIONAL CONSTRAINTS 

We assume a flat homogeneous and isotropic FRW universe such that 

= diag(l, a(t), a(t), a(t)), (7) 
where a(t) is the cosmological scale factor. By using 0, @, Q and Q we obtain 

T = -QH 2 , (8) 



where H = ^, being the Hubble parameter. The substitution of the vierbein (J7J) in (|6|) for 
i = = v yields 

T(l + 2f'(T))-f(T) = -lQnGp. (9) 
Besides, the equation i — 1 = v is 



|V=fflT/ w (r) + f'(T)[2T+ | J J>] - /(T) 



+ ^J^+T = lQnGp. (10) 



In equations rt9j) and (10), p and p are the dark matter energy density and pressure, respec- 



tively. It can be easily derived that they accomplish the conservation equation 



/— 3T 

p + y^-(p+p) = o. (ii) 

Here we assume that there are both matter and radiation components in the Universe, thus 
P = Pm + pr and p = |p r . If we rewrite the modified Friedmann equation (|9| in the standard 
form as that in general relativity, we can define a torsion energy density, expressed as 

PT= T ^(2Tf'(T)-f(T)). (12) 



An assumption that pr satisfies the torsion EoS, px = cotPt and by using (10) to define the 



torsion pressure pt, one can readily obtains the EoS parameter correspond to the DT as 

wr = -l-^, (13) 



where 



and 



M = (f(T) + 2Tf"(T))(l + 2f'(T) -121 + ?L) (14) 
B 1 = (1 + f'(T) + 2Tf"(T))(^p- - 2f(T)). (15) 



It is clear from ( 13 ) that the necessary condition for crossing the phantom divide line (ut = 
—1) is either 

f'(T)+2Tf"(T) = 0, (16) 

or 

n r = -3(l + 2/'(T)-^). (17) 

In addition, a dynamical EoS parameter also needs to satisfy the condition: 

d 



(U (Pt+Pt)^0, (18) 



when ut — > —1- Thus, one requires that the second term in the equation (13) changes sign 
in order for the parameter to exhibit a dynamical behavior. 

In the following, we are interested to examine the accelerated expansion universe and 
phantom crossing driven by torsion. For this we try with two /(T) models. We compare 
our both models with the ACDM model and also the SNe la observational data, combined 
with the information coming from the BAO (Baryon Acoustic Oscillation) and the CMB 
shift parameter. 

First Model: 

We consider the f(T) function as 

/(T) = a(-T)"sinh(^), (19) 
where a and n are two model parameters. This model is different from the one introduced in 



by multiplying f(T) by sinh(^) instead of tanh(^f). Substituting (19) into the modified 
Friedmann equation (|9]), we have 



a ~ (6# 2 ) n_1 [(2n - 1) sinh(l) - 2 cosh(l)] ' 1 ' 

Here Q m0 and Q r0 are the present dimensionless density parameters of matter and radiation, 



respectively. From f(T) given in (19), one can rewrite the DT energy density for this model 

as 

p T = a (~ T } n \(2n - 1) sinh(^) - 2^ cosh(^)]. (21) 



To satisfy energy condition, we must have px > 0. However, from numerical computation 
we obtain that for n > | the function < in some region. Therefore, in order to fulfill 
the energy condition requirement, the first constraint for n is n < 3 
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Also, to satisfy the equations (16) or (17), the necessary condition to cross the phantom 
divide line, —1, is 

tanh CBi] - {4n ~ 3)r ° r ■ (22) 



or 



Q r = -3 + 3a(-T) n - 2 ((l - 2n)Tsinh 

+2T cosh^J). (23) 

A dynamical dark torsion EoS parameter require to satisfy one of the above conditions in 
addition to the constraint (16) in our model. The given f(T) model satisfies the condition 
f(T)/T — > and therefore is consistent with the primordial nucleosynthesis and cosmic 
microwave background constraints [61]. Furthermore, for a — > 0, one regains the usual 
general relativity for an spatially flat FRW cosmology. 

Second Model: 

We consider f(T) as 

/(T) = a(-T)»coah(^)(l -exp(^)), (24) 

with three model parameters a, n and p. This model is also different from the one introduced 
in [60] by multiplying f(T) by cosh(^). By using (9), we can rewrite the parameter a in 
terms of other cosmological parameters as 

1 — ^m0 ~~ ^r0 



where 



(6# y^Ci v ; 



d = 2cosh(p)[pe p + n(l - e p )] + (e p - 1) x 

(cosh(p) + 2psinh(p)) . (26) 



From f(T) given in (24), one can rewrite the DT energy density for this model as 



Pt = [2cosh(— )[— e +n(l-e t )] 

+ (e t - i)[cosh( — ) -2— smh( — )]]. (27) 



Similar to the first model, for energy condition, px > 0, we need the constraint n < 1.44 



rounding up to 2 decimal. In order for ut crosses the phantom line, —1, the conditions (16) 



and (17) for this model become 



where 



and 



and 



tanh^V^f^ (28) 



A 2 = (nT 2 (l - 2n) - pT T(3 - An) - 4p 2 T 2 ) 

-(nT 2 (l - 2n) - 2p 2 T 2 ) (29) 



B, pTuiiin -3) + ( . X p( ^) (pT T(3-4n)+4p 2 T 2 ) (30) 



fi r = -3 + 6apT (-T) n ^ 2 1 - exp 



+3a(2n-l)(-Tr- 1 cosh(^ (l - exp . (31) 

Again, a dynamical dark torsion EoS parameter require to satisfy one of the above condi- 
tions in addition to constraint (16) for this model. Similar to the first model, the given 
f(T) satisfies the condition f(T)/T — > and therefore is consistent with the primordial 
nucleosynthesis and cosmic microwave background constraints. Note also that for a — > or 
p — > 0, the usual general relativity is retrieved. 



4. COSMOLOGICAL CONSTRAINTS AND TESTS 

In the following we first best fit the model parameters with observational data for lu- 
minosity distance of SNe la , the baryonic acoustic oscillation (BAO) distance ratio and 
the cosmic microwave background (CMB) radiation, using ?2 method. We then analyse the 
behavior of dark torsion EoS parameter and total EoS parameter predicted by model for 
the best fitted parameters. 
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4.1. Cosmological model constraints 



From observations, the Luminosity distance quantity, Dl(z), determines DE density. In 
the theoretical model with the H(z) obtained from numerical computation, we obtain 



D '«"< 1+ *>f^j- (32) 

where E(z) = H(z)/H . The difference between the absolute and apparent luminosity of 
a distance object is called distance modulus, fi(z), and given by, ji{z) = 5\og 10 D L (z) — /i 
where /i = 5log 10 h + 42.38 and h = (# /100)km/s/Mpc. 

Now, we obtain the constraints on the model parameters utilizing recent observational 
data, including SNe la which consists of 557 data points and belongs to the Union sample 
[75] , BAO distance ratio and CMB radiation. To constrain the parameters in the models 
from the SNe la, to best fit the observational data, one employs the x 2 value 



2 V- [/C^-^ri^r 

ASne — 2_^i ^2 ' y 66 ) 

i=l i 

where summation is over the cosmological data points. In (33), /i*' 16 and fi° bs are the distance 

modulus obtained from model and observation, respectively. Also, <jj is the estimated error 

of the fi° bs where obtained from observation. 

From joint analysis of the 2dF Galaxy Redshift Survey and SDSS data [T7] [78] , by using 

BAO data, the BAO distance ratio at z = 0.20 and z = 0.35 is obtained 

Dyiz = 0.35) . , 

-=P, { = 1.736 ± 0.065, 34 

D v {z = 0.20) ' v ' 

is a model independent measure with Dy{z) given by 

zbao , f ZBAO dz v 2i 1/3 



D ^^Hm { l m )] - (35) 

Therefore, one can obtain the constraint from BAO by performing the following \ 2 statistics 

[(D v (z = 0.35)/D v (z = 0.20)) - 1.736] 2 
Xbao ~ ^2 W 

We finally use the CMB data to constrain our model parameters. From CMB data, the 
CMB shift parameter R [79] [50]. contains the major observational information. Thus we use 
it to constrain the model parameters by minimizing 

2 _ [-R — Robs] 2 m \ 

ACMB — 2 ' \°') 

a R 




FIG. 1: The constraint on f2 m o and n at the 68.3%, 95.4% and 99.7% confidence levels from Sne 
la + BAO + CMB for the first model. 

where R b s = 1.725 ±0.018 [ST], is given by the WMAP7 data. Its corresponding theoretical 
value is defined as 

1/2 r ZcMB dz' 

R = Qmo J wr (38) 

with zcmb = 1091.3. The constraints from a combination of Sne la, BAO and CMB can be 
obtained by minimizing x| ne + Xbao + X 2 C mb- 

For the first model, from numerical computation, we find that the best fit values occur 
at Q m Q = 0.316 and n = 1.14 with Xmin = 560.86061. The contour diagrams at the 68.3%, 
95.4% and 99.7% confidence levels are given in FIG. (JTJ). From the graph, one conclude that 
with 68.3%, 95.4% and 99.7% confidence level the true values for both n and f2 m0 lie within 
the red, green and blue contours, respectively. For the second model, the best fit value of 
model parameters are fi m0 = 0.269, p = 0.05 and n = 1.12 with x 2 min = 544.20749. FIG. g 
shows the constraints on the model parameters at the 68.3%, 95.4% and 99.7% confidence 
levels . 

Alternatively, we can plot the likelihood for the pair model parameters in both cases(FIGS 
(|3j) and Q). Obviously in the first model since we only have one pair of parameters, there 
is only one likelihood plot for it. In the second model for the three parameters we have 
likelihood for three pairs of parameters. In FIG. Q, the bottom graph (the likelihood for 
the parameters p and f2 m o), one observes that there are two peaks, with the best fitted 
parameters are in the highest peak. 
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FIG. 2: The constraint on Q m o, n and p at the 68.3%, 95.4% and 99.7% confidence levels from Sne 
la + BAO + CMB for the second model. From top to bottom: (Q m Q-n) plane, (p-n) plane and 
(p-£l m o) plane. 




FIG. 3: Likelihood for the model parameters in the first model. 



The distance modulus, /i(-z), plotted in FIG. (|5]), for both models with the best fitted 
parameter values are compared to the ACDM model by combining Sne la, BAO and CMB 
observational data. 

4.2. Dark torsion EoS parameter 

Using the best-fitted constrained model parameters in the last section, we can now 
discuss the behavior of dark torsion EoS parameters in the models. In the first model, FIG. 



II 




FIG. 4: Likelihood for the model parameters in the second model. 




FIG. 5: The distance modulus, fi(z), with respect to z expected for the standard concordance 
model ACDM with = 0.74 and O m o = 0.26 (solid line), and for our models (dashed line) with 
the values of the best-fit coming from Snla+BAO+CMB. Top panel: the first model. Bottom 
panel: the second model. 

(JgJ) , bottom), shows the dynamics of the dimensionless density parameters for radiation, 
matter and DT in the universe. It reveals that the universe in this model has a long 
enough period of radiation domination to give the correct primordial nucleosynthesis and 
radiation-matter equality and a matter dominated phase. In other words, the usual early 
universe behavior can be successfully obtained to agree with the primordial nucleosynthesis 
and the cosmic microwave background constraints. It also shows that for this model the 
DT start dominating the universe at about z ~ 0.23. 

In addition, in FIG. ^ top), the evolutionary curves of the dark torsion EoS parameter 
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FIG. 6: The trajectories of the dark torsion EoS parameter (top panel) and the trajectories of the 
dimensionless density parameters for the DT (green), radiation (red) and matter (black) (bottom 
panel) with the best fitted values of n and f2 m o in the fist model. 

for the best fitted values of n and Q m Q are shown. One can see that the dark torsion EoS 
parameter crosses the phantom divide line from the values greater than —1 (non-phantom 
phase) to smaller than —1 (phantom phase) in the past and become tangent to —1 in the 
future. The result is in consistent with the one obtained in [64J where the universe in the 
f(T) theory finally enters a de Sitter expansion phase. Also, it can be easily shown that the 



crossing occurs at the red shift z where Eq.(22) is satisfied and the constraint (23) is related 
to entering a de Sitter expansion phase. 

For the second model, in FIG. ([7]) bottom) we show the trajectories of the dimensionless 
density parameters and the dark torsion EoS parameter for the best fitted model parameters 
of p, n, and Q m o and the DT start dominating the universe at about z ~ 0.4. One can see in 
the top panel that the dark torsion EoS parameter crosses the phantom divide line from the 
values greater than —1 to smaller than —1 in the future and finally become tangent to —1 in 



the near future. Also, in this figure one can check that the conditions (30) or (31) satisfies 
when the crossing occurs in the early future and the model enters a de Sitter expansion 
phase in the late future. 
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FIG. 7: The trajectories of the dark torsion EoS parameter (top panel) and the trajectories of the 
dimensionless density parameters for the DT (green), radiation (red) and matter (black) (bottom 
panel) with the best fitted values of n, £l m o and p in the second model. 

4.3. Total EoS parameter 

In the pervious section we examined the dynamics of the dark torsion EoS parameter 
which in the first model shows phantom crossing in the past plus de Sitter phase in the 
future and in the second model shows phantom crossing in the early future and de Sitter 
phase in the late future. However, the universe is filled with the CDM and radiation and a 
physical model requires the inclusion of correction terms in the formulation. One can rewrite 
the Friedmann equations (|9| and (10) as 



T = -16nGp tot (39) 

2 7=if + T = 167rGKt (40) 

where p to t = p + Pt and p to t = P + Pt and the energy density, p, is for CDM and radiation 
filled the universe. Simply, the total EoS parameter in terms of the redshift z as 

1 (l + z)dT 

Using the best fitted model parameters obtained from \ 2 method, one observes the evolution 
of the total EoS parameter bj tot as a function of z for the above two models. 

FIG. (Jsj) top) shows the last three phases of the evolution of the universe, i.e., the radia- 
tion dominated (u to t = 1/3), the matter dominated {u to t = 0) and the late time acceleration 
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FIG. 8: The trajectories correspond to the total EoS parameter as a function of z expected for the 
standard concordance model ACDM with £7 a = 0.74 and £l m o = 0.26 (dashed line), and for our 
models (solid lines) with the values of the best-fit coming from SNIa+BAO+CMB, for the first 
model in two different ranges 

for the first model in comparison with the ACDM model. In more details to observe the 
behavior of the u to t m the near past and future, the FIG. (|8]) bottom) shows that the uni- 
verse transits from deceleration era to acceleration era when u to t — — 1/3 at about z ~ 0.61 
and approaches a de Sitter phase in the future where Utot — 1- We also see that the cur- 
rent value of u^t at z = is about Utot — —0.82 whereas for ACDM model is about ~ —0.74. 

Also, FIG. Q top) shows the last three phases of the evolution of the universe, i.e., 
the radiation dominated (u) to t = 1/3), the matter dominated (u tot = 0) and the late time 
acceleration for the second model in comparison with the ACDM model. For the behavior 
of the Utot in the near past and future, the FIG. (J9j) bottom) shows that the universe 
transits from deceleration era to acceleration era when utot = —1/3 at about z ~ 0.74 and 
approaches a de Sitter phase in the future where u to t — > — 1- The current value of u to t in 
this model is about u to t — —0.73 whereas for ACDM model it is about ~ —0.74. 



5. SUMMARY AND REMARKS 

In this manuscript, the theory of T+f(T) based on TEGR where torsion is the geometric 
object describing gravity instead of the curvature and the 2nd order differential equations 
that is remarkably simpler than f(R) theories is investigated for two different f(T) models. 

15 




FIG. 9: The trajectories correspond to the total EoS parameter as a function of z expected for the 
standard concordance model ACDM with Q\ = 0.74 and £l m o = 0.26 (dashed line), and for our 
models (solid lines) with the values of the best-fit coming from SNIa+BAO+CMB, for the second 
model in two different ranges. 

Our purpose in this work is, by using numerical methods, to reproduce the recently detected 
acceleration of the universe and the dynamics of the EoS parameter for DT without resorting 
to the DE. Our work is distinguished from the work given in [60] by an analytic discussion 
on model equations for phantom crossing behavior for the given f(T) models. In addition 
the parameters in the model are best fitted with the observational data in advance, before 
numerical computation of the EoS parameters for the models. This gives us a more accurate 
picture of the universe and physically meaningful solution for the derived parameters. 

Our original motivation for studying f(T) gravity and considering the specific models as 
an alternative to conventional dark energy in general relativistic cosmology is to interpret the 
current universe acceleration and reproduce a dynamical dark energy in agreement with the 
observations. Our f(T) models satisfy the conditions at the high redshift, consistent with 
the primordial nucleosynthesis and cosmic microwave background constraints. Furthermore, 
the usual general relativity can be retrieved at low energy limit. 

We tested the observational viability of the models by using the most recent SNe la data 
combined with the information coming from BAO peak and CMB shift parameter in order 
to find constraints on the model parameters, n, p and fi m0 . The numerical analysis indicates 
that these two types of f(T) theories are compatible with observations. At 68.3% confidence 
for the first model we found that the parameter's values lie in the ranges n = [0.9, 1.3] 
and fl m o = [0.288,0.345]. The model with the best fitted values minimizing the \ 2 that 
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combines SNe Ia+BAO+CMB data given by n = 1.14 and Q m o = 0.316. Furthermore, 
for the second model, at 68.3% confidence we found that the parameter's values lie in the 
ranges n = [0.79,1.44], Q m0 = [0.243,0.302] and p = [-0.48,0.17]. The model with the 
best-fit values minimizing the x 2 that combines SNe Ia+BAO+CMB data for this model 
given by n = 1.12, p = 0.05 and f2 m0 = 0.269. 

With the best fitted model parameters for the distance modulus and in comparison with 
the observational data we calculate the dark torsion EoS parameter and also dimensionless 
density parameters for the two models. In the first model, the phantom crossing occurs in 
the past, while currently in phantom era and approaches a de Sitter phase in future when 
z — — 1 and universe become infinite. For the second model the universe is currently in 
quintessence era, undergoes phantom crossing in near future and again approaches a de 
Sitter phase later in future when z — — 1 and universe become infinite. 

However, in reality, in addition to the repulsive DT, the universe is also filled with CDM 
and radiation. Thus the dynamics of the universe has to be explained with the total EoS 
parameter, u tot . The computation given for u tot exhibits the last three phases of cosmo logical 
evolution: radiation era, matter era and late acceleration where in the last stage transition 
from deceleration to acceleration occurs for the first model at z ~ 0.61, the second model 
at z ~ 0.74 and for ACDM model at z ~ 0.75. Also, from the graphs the current values 
for u to t for the first model is uj to t ~ —0.82 and the second model is u tot ~ —0.73 whereas 
for the ACDM model is ~ —0.74. One can conclude from findings that the second model 
is more compatible with the ACDM model [H2]- We also find that the values obtained for 
the total EoS parameters in both models are within the la confidence region of current EoS 
parameter in the WMAP 5-year results [83] . 
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